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^ . Abstract 

The one-dimensional (ID) t-J model is investigated by using a Gutzwiller-Jastrow-type 
variation method and the exact diagonalization of small systems. Variational expecta- 
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tion values are estimated by the variational Monte Carlo method with sufficient accu- 
racy. First, we represent the diagonalization results. Physical quantities like momentum 
distribution and some correlation functions show some behaviors which are not expected 
in repulsive models, as the value of J/t increases. These properties as well as energy are 
well understood by introducing intersite correlation factors into wave functions. The 
phase transition to a separated phase in large- J/t region can be described by an attrac- 
tive Jastrow wave function in quantitative agreement with the exact results. On the 
other hand, for the supersymmetric case (J/t = 2) the original Gutzwiller wave function 
becomes an extremely good trial function for all the range of electron density. Here a 
kind of "free electron" state is realized, particularly in the low electron density. Next, 
the above wave functions are compared with the Tomonaga-Luttinger-liquid-type wave 
function proposed by Hellberg and Mele. It is found that the correlation factors in short 
distances control bulk quantities like energy and the magnitude of correlation functions, 
while the long-range part of correlataion factors determines the critical behavior of cor- 
relation functions. Lastly, using these functions, charge and spin susceptibilities and 
magnetization curve are estimated, which agree with the exact results. It is shown 
that the Mott transition in ID t-J model is quite different from the Brinkman-Rice 
transition. 

PACS numbers: 75.10.Jm, 71.10,+x, 74.65. +n, 74.70.Vy 
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§1. Introduction 

The t-J model is an important model to study highly correlated electron systems for 
its simplicity and close relationship to the high temperature superconductivity. 1 -* Many 
properties in one-dimensional (ID) systems have been clarified extensively by a number 
of methods: Bethe-Ansatz solutions, g-ology, Tomonaga-Luttinger (TL) liquid theory, 
quantum Monte Carlo simulations, exact diagonalization studies of small clusters, and 
conformal field theory. We expect that the study of ID systems will shed light on more 
realistic higher-dimensional systems and that the comparison of the various methods 
with the well-established ID results enables us to judge the validity of such methods 
and approximations. 

In contrast to the ID Hubbard model, 2 -* the Bethe Ansatz solution does not exist 
in the ID t-J model except for J/t = (spinless-fermion case) and J/t = 2 (super- 
symmetric case). 3 ' 4 ) In both soluble cases, TL liquid 5 '' is realized and the exponents of 
long-range behaviors of correlation functions were calculated exactly by combining the 
Bethe Ansatz equations and the conformal field theory. 6-9 ^ Also obtained were bulk 
quantities like spin susceptibility % s , charge susceptibility % c , specific heat coefficient 
and effective transport mass, which characterize metal-insulator (Mott) transitions. 10 -* 
On the other hand, for general values of J/t, Ogata et al. 11 -* studied the low- lying energy 
spectrum of finite systems to obtain the correlation exponents. In the phase diagram of 
J/t and the electron density n, the TL liquid theory holds in the small- J region below 
J c /t = 2.5 ~ 3.5, depending on n. A phase separation takes place in the larger-J/t; 
there is a region in which the superconducting correlation is dominant, between the 
phase separation and the supersymmetric case. 

Meanwhile, as for physical quantities such as momentum distribution function or 
spin- and charge-correlation functions, only the long-range behaviors were clarified in 
the above analytic methods; the global features were calculated numerically. In the 
limit of J/t — > 0, identical with the \arge-U limit of the Hubbard model, the correlation 
functions were obtained by taking advantage of the spin-charge separation in the ground 
state. 12 -* For the other values of J/t, Assaad and Wiirtz 13 ) and Imada 14 -* have carried out 
quantum Monte Carlo simulations. Pruschke and Shiba 15 -* studied the superconducting 
correlation functions by the exact diagonalization. All these results are consistent with 
the correlation exponents obtained by the analytic methods. 

Although the ground-state properties in the ID t-J model have been clarified quite 
well, it is still important to examine variational wave functions. For the explicit form 
of the wave function will make the complicated physics easy to grasp. So far, various 
kinds of variational states have been proposed for strongly correlated systems. 16_23 ) The 
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Gutzwiller wave function (GWF) 24 ) was studied numerically 16 ' 17 ) and analytically. 18 ' 19 ) 
These studies concluded that the GWF is excellent for the one-dimensional Heisenberg 
model, but is unsatisfactory, even qualitatively, in describing the properties of the strong 
coupling Hubbard model or of the small-J region of the t-J model. For example, the 
GWF does not reproduce the 2kp peak in the spin correlation function; in the momen- 
tum distribution, it has a strange enhancement for k > k-p. The main reason is that 
the density correlation is not sufficiently introduced in the GWF, although the spin 
correlation is well incorporated. These unsatisfactory features are partly remedied by 
introducing Jastrow-type intersite correlation factors. 20 ' 22 ' 23 ) 

In this paper we show that the ground state properties of the ID t-J model become 
more easily accessible from a variational viewpoint. First, the behaviors obtained by 
the exact diagonalization for J/t = 2 is described extremely well by the GWF. Next, 
the wave function is improved for other values of J/t by introducing Jastrow-type wave 
functions with intersite correlation factors, which are classified into two types: one is the 
conventional repulsive and attractive Jastrow-correlation factors and the other is a long- 
range Jastrow factor introduced by Hellberg and Mele. 23 ) The exponent in the former 
state is the same as the Fermi liquid, while the latter state has non-trivial exponents. 
It is found that the long-range behavior of the Jastrow factor is essential for the non- 
trivial exponent, or non-Fermi liquid behavior. On the other hand, the variational 
energy and the global features of the correlation functions are modified mainly by the 
short-range behavior of the Jastrow factor. We also show that quantities like Xc Xs and 
magnetization curve obtained by the above wave functions are not only qualitatively 
but quantitatively consistent with the exact results. This aspect is in sharp contrast 
with the so-called Brinkman-Rice transition. 

The outline of this paper is as follows: in §2 we discuss the diagonalization results 
for various correlation functions as well as the ground-state energy. In §3 the prop- 
erties of the GWF are studied comparing with those of the exact calculations for the 
supersymmetric case. Section 4 is devoted to the study of the other values of J/t by 
introducing Fermi-liquid-type Jastrow wave functions. In §5 another variational state 
with an essentially long-range Jastrow correlation (TL-liquid-type) is examined and 
compared with the former type. In §6, Xc Xs and magnetization curve are investigated. 
Section 7 is assigned to summary. In Appendices A and B, an analytical approach to 
the GWF used in §3 and the behavior of Xc and Xs in the Gutzwiller approximation 
compared in §6 are summarized respectively. A part of the results in this paper was 
reported before. 22 > 25 ) 
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§2. Ground-State Properties Obtained in Small Clusters 

We study the one-dimensional (ID) t-J model denned as 

H = H t + Hj, (2.1) 
H t = -t J2(c} a c j+la + h.c.), (2.1a) 

jo 

Hj = J J2( S j ' S J+ 1 ~ ^ n J n i+i)' ( 2 - lb ) 
j 

in the subspace with no double occupancy with t, J > 0. Spin operators vanish when 
they are applied to empty sites. Henceforth we take t as the unit of energy. 

We use the Lanczos method and the conjugate gradient method 26 ^ to obtain the 
ground-state wave function in small clusters. Figures l(a)-(c) show the momentum 
distribution function 

n(k) = (c{ a c ka ), (2.2) 
and spin- and charge-correlation functions 

^) = ^rE 4 ^^) elfc(rj " r ' ) ' 
iVa j,i 

N ^ = W E^K'^) " (nj)(n e )}e ik ^-^, (2.3) 



obtained exactly for various values of J. These results are for the quarter-filled case: n = 
N/N a = 0.5, with N and N a being the number of electrons and sites, respectively. As 
shown in the ID U — > oo Hubbard model, 12 ^ even the 16-site cluster gives fairly reliable 
correlation functions. The thermodynamic limit is estimated easily by interpolating the 
data points for the 16-site cluster except for the singularities. Notice that S(k) and N(k) 
are also calculated by Assaad and Wiirtz 13 ) using quantum Monte Carlo simulations. 
Their results agree with Figs. l(a)-(c). 

There are power-law singularities in the momentum distribution n(k) at /cf and 
3/cf- The 3/cf singularity is not so prominent in the small systems in Fig. 1(a). The 
power-law singularity at k-p is the strongest when J/t = 0. With increasing J/t, its 
singularity becomes weaker until J/t = 2 where n(k) seems to have a jump at k-p. 
In the region J/t > 2 the power-law singularity shows up again. This dependence 
on J is consistent with the behavior of the numerically obtained correlation exponent 
K p which increases from 1/2 as J increases. n ) Here the exponent K p is defined such 
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that (sfsj) decays as exp{2ik F rij} /r^ p . The exponent at k F of n(k) is defined as 
n(k) = n — Csgn(/c — k F )\k — k F \ a . Since this exponent a is related to K p as a = 
(Kp + V-Kp — 2)/4, it starts from | at J/t = 0, decreases for < J/t < 2, becomes 
a = near J/t = 2 and then begins to increase for J/t > 2 where K p becomes larger 
than 1. 

There are 2k F singularities in S(k) and N(k). The TL liquid theory predicts that 
these singularities have the same correlation exponents K p and that the size-dependence 
is A-CN a As seen from Fig. 1, it is difficult to estimate K p from the size-dependence 
of the 2k F peak. Especially the peak in N(k) is much weaker than that in S(k). We 
expect the coefficient C of the size-dependence N a Kp is very small for N(k). 

The global features of the correlation functions for small values of J resemble the 
results in the large-t/ Hubbard model. 12 ) As J increases they lose this behavior and near 
the supersymmetric case (J/t = 2), the system behaves similarly to the non-interacting 
case. This again corresponds to the fact that the exponent K p becomes f (free-electron 
value) near J/t = 2 (J/t ~ 2.3 for n = 0.5). 11 ) N(k) and S(k) become almost flat in the 
region k > 2k F , which is the same behavior as in the non-interacting case. However, 
note that this global features of correlation functions are non-trivial even if K p = f, 
since this exponent only guarantees that the long-range behavior of correlation functions 
is the same as in the non-interacting case. Actually the absolute value in the flat region 
is quite different from the non-interacting value. For the case of 8 electrons in f 6 sites, 
we get 

N(k > 2k F ) = 0.312 ~ 0.3f9, 

S(k > 2k F ) = 0.698, (2.4) 

while in the non-interacting case, N(k > 2kp) = S(k > 2/c F ) = n = 0.5. The sum of 
these two values is, however, 

N(k > 2k F ) + S(k > 2k F ) = 1.001 ~ 1.009, (2.5) 

and this is surprisingly close to the non-interacting value (2n = 1). We will discuss this 
result in the light of the GWF in the following chapter. 

Next we estimate the ground-state energy. The ground state is always singlet and 
non-degenerate, if we choose periodic (antiperiodic) boundary conditions for N/2 = 
odd (even), respectively. 12 ' 11 ) Under these boundary conditions, the energy converges 
smoothly to the thermodynamic limit. For n = 0.5 we calculate the ground-state 
energies in 4, 8, 12, and 16 sites clusters and fit the results to the formula 

E/N a = Coo + C^Nl + C 2 /Ni + C 3 /iV a 6 . (2.6) 
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In Fig. 2 the fitted values of are shown by a solid line. In the region J/t > 3.4 
the energies cannot be fitted to this formula because the system phase separates in this 
region and the size-dependence is different from (2.6). To check the convergence for 
J/t < 3.4, we calculate another series of singlet energies by using a different boundary 
conditions, i.e., antiperiodic ones for N/2 = odd and vice versa. 27 ) Fitting of the data to 
the same formula (2.6) gives another estimate of the ground-state energy, (e^J. [In this 
case Ci > 0, while the former fitting gives C\ < 0.] The fitted values are shown in Table 
I. The difference between and e'^ is very small (A eoo < 10" 6 t), 28 ) which ensures the 
reliability of the estimate. For n = 3/4, we estimate similarly the ground-state energy 
by using 8 and 16 sites clusters with 6 and 12 electrons, respectively. The obtained 

is shown in Table II. Because in this case we fit the data with two parameters 
and Ci, the error is larger than that in the quarter-filled case. We will use later to 
compare with the variational energies. 

§3. Comparison with the GWF Near the Supersymmetric Case 

In this section we compare the results in the previous section with the Gutzwiller wave 
function (GWF). The GWF is defined as 

PA=n(i-wi^, (3.i) 

3 

where $f is a simple Fermi sea. In ID case, the analytic expression for the physical 
quantities were developed. 18 > 19 ) 

First let us compare the variational energy. The expectation value of the kinetic 
energy is calculated as 



Et = ^ = fjk cos k(clc k(7 ), (3.2) 



where (• • •) indicates the expectation value in the GWF. The analytic expression for 
{ c \ a c ka) has been given by an infinite summation. 18 ) The detailed calculation is sum- 
marized in Appendix A. Here we show that the exchange energy can be obtained in a 
compact form. 19 ) From the expression of S(k) and N(k) (see below and Appendix A), 
we get 

(4^ +1 ) = -i{Si(7r)-Si((l-n)7r)}, 
I 

(nin i+1 ) =n 2 + ( cos 2mr ~ 1) (3.3) 



+ - (^HHL + (1 - n ) cosmr^) {Si(yr) - Si((l - ra)*-)}, 

71 \ 71 J 
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with n = N/N a being the electron density. Using a fact that the GWF is singlet, 29 ) we 
obtain analytically 



JVa 

4tt\ 



si n 11 it 

{3 + -—- + (1 - n) cosn7r}{Si(7r) - Si((l - n)n)} 



+ -nn 2 + -!-(cos2mr - 1) 1. (3.4) 

Z7T J 

In Fig. 3, we compare the total energy of the GWF with that of the Bethe Ansatz 
(BA) at J/t = 2. It is surprising that two results are indistinguishable for any value of 
n in the scale of this figure. In fact, for the half-filled or the Heisenberg case (n = 1), 
the energy of the GWF, 16 ' 18 ) 

E(GWF)/t = - i[— Si(yr) + 1] = -1.384235 • • • , (3.5) 

2 TV 

is extremely close to that of the Bethe-Ansatz result, 

£7(13 A) ft = -2 In 2 = -1.386294 • • • . (3.6) 
For n = 0.5, analytic expressions (3.2) and (3.4) give 

E t /t = -0.574632 ••• , (3.7) 



Ej/J= -0.164230 (3.8) 
Thus the total energy of the GWF becomes 

E(GWF)/t= -0.903092 (3.9) 

which is quite close to the exact one: 

E{BA)/t= -0.903649 (3.10) 

The difference is only 0.06%, which is better than the half-filled case (0.15%). 

In the low density region, we find that the coincidence is much better. We can see 

that the ground-state wave function for 2-electron problem at J/t = 2 is given by 



l*> = M=o T 4=04lO>- 
7 



(3.11) 



This means that the GWF is exact in the low density limit. In fact, the total energies 
of the GWF at J/t = 2 is 

7T 2 

E(GWF)/t = -2n+ —n 3 +0(n 4 ), (3.12) 

J. — ■ 

which coincides with the BA results up to the order of n 3 . The detailed calculation 
for (3.12) is also given in Appendix A. This is consistent with the fact that the critical 
exponent K p approaches 1 for n — > 0. This Fermi liquid state is nothing but the GWF. 

In Fig. 4, we compare n(k) for the GWF with that for the small cluster calculation. 
As for n(k), two results almost coincide quantitatively except for the singularity at k-p. 
It is remarkable that the ground state of the t-J model has an enhancement of n(k) in 
the vicinity of 7r, which was considered before as a pathological behavior of the GWF. 17 ' 
It can be shown that this enhancement originates from the correlated electron motion, 

(4* c jM ~ ni-*)(l - %-a))o, (3.13) 

where (• • -)o represents the expectation value in the free Fermi see $p without the 
Gutzwiller projection. After a straightforward calculation, (3.13) becomes 

PijaO- ~ n-a + n 2 _ a + 5 lJ n- a - Pij^Pji^), (3.14) 
and its Fourier transform is 

n a (k)-2n_ a n a (k) + f 2 (k). (3.15) 

Here Pij a is the Fourier transform of n° a (k) = 6(kp — \k\), and f 2 is given in Appendix A 
as one of the lowest order terms in the analytic calculation. The enhancement of n(k) 
in the region of k > kp is roughly reproduced from this simple lowest-order, (3.13). 

The slight difference between the GWF and the exact ground state is the behavior 
around kp which originates from the TL liquid nature. The GWF is essentially a Fermi 
liquid state with discontinuity q = ^/l — n at k = kp. On the other hand, n(k) in the 
true ground state has a power-law singularity as 

n(k) = n(kp) — c\k — kF\ a sgn(k — fcp) (3.16) 

with a = (K p + 1/Kp — 2)/4. For n = 0.5, K p is about 0.85 and is not far from 1 of the 
Fermi liquid case. 11 ' 13 -' This is the reason for the two results are relatively close to each 
other even in the vicinity of k-p. In order to describe these power-law behaviors around 
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k F , it seems necessary to introduce low energy excited states around the Fermi surface 
into the trial wave function, 20 ) or a long-range Jastrow factor 21 ' 23 ) as discussed in §5. 

In Fig. 5, S(k) and N(k) are compared with the GWF. S(k) increases for k < 2k F 
and becomes constant for k > 2k F : 



I -ln(l -n) , 
(2k F = nir). N(k) has different curves 



< k < 2k F 
2k F < k < n 



(3.17) 



i \k\ ln 



7T 



(1-n) ' 

(0 < k < min(2/c F , 2tt - 4fc F )) 
^rf) + ln(l + n-M), 

(1— n) \ 7r / ' 

(2A; F < < min(4/c F , 2n - 4k F ) 
when 2/c F < 2n - 4k F , or n < 2/3) 



2 - 2n + 



\k\ 



ln 



iV(fc) = < 



2- ^ (l-n+M) ' 

(2n - 4k F < k < 2k F 



when 2tc - 4k F < 2k F , or 2/3 < n < 1) 



(3.18) 



2n + ln(l - n) , 



(Ak F < k < tv 

when 4/c F < 7r,or n < 1/2) 



2 - 2n + ^ ln 7 



Mi.(!^±l) 



l+n- 



+ ln(l + n-M) 



(max(2/c F , 27r — Ak F ) < k < n 
when 7r < 4k F , or 1/2 < n < 1) 



Figures 4 and 5 show that the GWF is an extremely good trial function to describe 
globally the ID t-J model in the supersymmetric case, although it is not the exact 
ground state. The J-dependence of the variational energy for the GWF is also shown 
in Fig. 2 using (3.7) and (3.8). It is apparent that the GWF is very good only in the 
vicinity of the supersymmetric case. 

It is possible to speculate the reason. An up-spin electron, for example, can hop 
to a neighboring site with the same magnitude t = J/2 independent of the state of 
the neighboring site (hole or down-spin electron). 4 ) This fact makes the up-spin part of 
the wave function similar to the Fermi sea where the up-spins move freely without any 
correlation to the down-spins. Actually the electron system without any interaction has 
an energy 



4 nn 
E t /t = sin — 

7T Z 



-2n + ^-n 3 +0(n 5 ), 
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which coincides with the exact result up to 0(n s ). 

Now we discuss the flat region of S(k) + N(k) in k > 2kp which we have shown in 
the previous chapter (Fig. 1). The summation of them is rewritten as 



Therefore it represents the density correlation between the same species of spin. The 
coincidence of this quantity to the non-interacting value shows that the each spin be- 
haves freely. The Gutzwiller projection only affects the correlation between the different 
species of spins, such as {n^n^). The fact that the GWF behaves almost as a free state 
will be confirmed also for such quantities as charge and spin susceptibilities and mag- 
netization curve in §6. 

It seems that the fact the GWF is excellent in the supersymmetric case does not 
occur by chance. Recall that for the ID supersymmetric t-J model with 1/r 2 long-range 
coupling, the exact ground state is the GWF. 30,31 ' The lowest eigenvalue specified by 
electron density n and magnetic polarization m = (N^ —Ni)/N a consists of independent 
contributions of n and m; the logarithmic correction does not exist for the critical 
behaviors. It is surprising that except for the long-range behaviors the main features of 
this long-range model are shared with the ordinary t-J model. Furthermore recently it 
has been found that many low- lying excited states of the above long-range t-J model, 
including the case of finite magnetic polarization, are represented also by Gutzwiller- 
projected determinantal wave functions. 32 ' 25 ' 33 ^ Thus to elucidate the relation between 
these models and to study excited states by the variational approach are fascinating 
future problems. 
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N a 



(3.19) 
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§4. Fermi-Liquid-Type Jastrow Functions 

In this section we study the Fermi-liquid-type wave functions to describe the t-J model 
away from J ft = 2, keeping the application to higher dimensions in mind. In many- 
body problems Jastrow-type wave functions with two-body correlation factors are often 
used. First, notice that the trial state becomes non-singlet if the Jastrow correlation is 
spin-dependent. 34 ) Therefore we study here only the spin-independent charge density 
correlation: 35 ^ 

* = II II i 1 - i 1 - V(r je ))n ja nea'] $f, (4.1) 

jt era' 

where rjg = \rj — rg\. For rj(r) simple forms are desirable, satisfying the condition of 
the t-J model: rj(0) = 0. This condition is to project out the double occupancies. In 
this section we consider two cases in addition to the GWF: 

r J arctan J (a) (RJWF) 

V(r) = I 1 ' (b) (GWF) , (4.2) 

1 1 + ^ (c) (AJWF) 

where r ^ and £, a and (3 are positive variational parameters. The typical forms of 
rj(r) are shown in Fig. 6. In these function, 77(00) = 1 holds and the value of rj(l)/rj(oo) 
is finite for every parameter set. This point is essentially different from the function 
with a long-range Jastrow factor discussed in the next section. 

The form (a), which we call the repulsive Jastrow wave function (RJWF) in this 
paper, includes repulsive correlation. It prefers configurations with electrons mutually 
apart. In the limit £ — > 0, it is reduced to the GWF (b). This wave function was 
previously introduced to study the repulsive Hubbard model in strong correlation. 20 - ) In 
that work it was found that the RJWF lowers the variational energy in the Hubbard 
model and it reproduces qualitatively the enhancement at 2/cp of S(k). Therefore we 
expect that the RJWF is suitable also for the small-J region of the t-J model. 

An attractive Jastrow wave function (AJWF) with the correlation factor (c) favors 
local configurations with electrons close to each other. The parameter a adjusts the 
amplitude of such attractiveness; as a increases, more emphasis are laid upon attractive 
electron configurations. On the other hand, the other parameter (3 controls the decaying 
behavior of 77(7"). With increasing /3, effective attractive range becomes narrower. It is 
reduced to the GWF when a — > or f3 — > 0. We will show that the AJWF properly 
describes a homogeneous phase for relatively small a as well as a fully phase-separated 
state for large values of a. 

Comparing with the RJWF and the AJWF, we may regard the GWF as a "free- 
electron" state in that there is no amplitude modification from the non-interacting state 
except for the exclusion of doubly occupied sites. 
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We evaluate the expectation values by the variational Monte Carlo (VMC) method. 16 ' 
This method gives virtually exact expectation values for given trial functions. In the 
VMC calculations in this work we use systems with electron number N = 41 + 2 (I: 
integer) with the periodic boundary condition. Sample numbers (3 x 10 4 ~2x 10 5 ) and 
sampling intervals (50 Monte Carlo steps at the maximum) are taken so as to reduce 
statistical fluctuations enough. In this section we use typically 60- and 72-site systems. 

Before going to the cases of J/t ^ 2, we check the supersymmetric case. Total 
energy per site of the RJWF and the AJWF for J/t = 2 and n = 0.5 are shown in Figs. 
7(a) and 7(b), respectively. These figures tell us that the Jastrow correlation factors 
raise the variational energy and that the GWF is the best trial state for J/t = 2. In 
the following we will see that the RJWF has lower energy than the GWF for J/t < 2 
(§4.1) and the AJWF has lower energy than the GWF for J/t > 2 (§4.2). 

4-1 Repulsive Case 

We starts with a brief survey of energy. In Fig. 8, Et/t and Ej/J in the RJWF are 
plotted as a function of ( for n = 0.5. As ( increases, electrons keep apart from each 
other and thus E t tends to decreases and Ej increases. We show the total energy 
E(Q for J/t = 1.0, as an example, in Fig. 9 using the values in Fig. 8. The minimum 
is situated at £ ~ 0.7; the energy is considerably improved upon the GWF (£ = 0). 
Similarly we find that the minimum appears at ( ~ 1.6 (for J/t = 0), 1.2 (0.5), 0.7 
(1.0), 0.4 (1.5), and (2.0), respectively. The repulsive interaction represented by the 
magnitude of £ becomes weaker with increasing J. Variational energies thus obtained 
are plotted also in Fig. 2 to compare with the exact ground-state energy. 

Next, let us consider the critical value J c below which the RJWF has a lower 
variational energy than the GWF. In order to estimate J c we only have to see the in- 
dependence of the variational energy near £ = 0: if the slope of the energy at £ = is 
positive, the optimal state is £ = 0, which is the GWF. On the other hand, if the slope 
at £ = is negative, the variational energy has a minimum at a certain value of £ > 
and the RJWF is more stable than the GWF (C = 0). 

In Fig. 10 variational expectation values, Et/t and Ej/ J, in the RJWF are shown 
for n = 0.3 and small £. From this figure we read Et(()/t = —0.471 — 0.0413^ and 
Ej(()/ J = —0.0540 + 0.0217C- Since the total variational energy is given by 

E(0=E t (0+Ej(C), (4.3) 

the slope of E(() near £ = is -0.0413* + 0.0217J. This means that for J/t < 
0.0413/0.0217 ~ 1.9, the slope is negative and the RJWF improves the variational 
energy upon the GWF. In a similar way, we obtain the value of J c /t also for n = 0.5, 
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0.75 and 0.833 • • •. As n increases, the statistical fluctuation becomes severe and it is 
not easy to estimate accurate values. However, we find J c /t = 2.0 ± 0.2 for the above 
densities. From this we conclude that a repulsive correlation is necessary to reproduce 
the ground state for J/t < 2. Actually the value of J c /t ought to be slightly shifted to 
a larger value of J/t in the high density region, as will be seen in the next section. This 
is consistent with the diagonalization study. 11 ) 

Leaving the discussion of correlation functions in the next section, here we focus on 
the discontinuity q of momentum distribution n(k) at /cf, which is one of the character- 
istic properties of a Fermi-liquid-type wave function. This quantity is equal to the wave 
function renomalization factor and related to the inverse of the effective mass m* /m in 
the Fermi liquid theory. Similarly it is the energy reduction factor in the Gutzwiller 
approximation. In Fig. 11 the value of q is plotted as a function of n calculated with 
the optimized wave function for each value of J/t. 

For the supersymmetric case the result of the GWF, 18 ' q = VT — n, is shown. In 
this case q becomes 1 in the limit of n — > 0, which is the value of the non-interacting 
system. In the other limit n — > 1, q vanishes. This implies the Fermi surface vanishes, 
namely effective mass diverges, corresponding to the metal-insulator transition. As J/t 
decreases, q becomes small; the correlation affects the whole electron density in this 
case. Although the result for J/t = 2.5 includes relatively large statistical error, it is 
obvious that q tends to decrease, with increasing J/t from J/t = 2.0. Thus we conclude 
that the supersymmetric case is the most weakly interacting; as J/t goes away from it, 
the correlation effect becomes severer. 

4-2 Attractive Case 

We study the AJWF in this subsection. For J/t > 2, it is natural to expect that the 
attractive correlation between electrons is dominant. In Fig. 12 we show the variational 
expectation values of Tit and Tij for the AJWF (n = 0.75). Ej becomes lower and Et 
becomes higher as a increases, because the amplitude of configurations with electrons 
located next to one another increases. Note that E t abruptly approaches to zero and 
Ej to the energy of the spin system in the region a > 10 for (3 = 0.625 ~ 1. This is 
nothing but a sign of the phase transition. We will come back to this point shortly. 

First we estimate the critical values above which the GWF becomes unstable against 
the AJWF by studying the small-a (and various f3) behavior, as in the RJWF case. We 
obtain J c /t = 2.0 ± 0.1 for n = 0.5 and 2.0 ± 0.4 for n = 0.25 and 0.75. 

Searching the energy minimum in the a-f3 plane, we obtain the optimized energy 
for J/t > 2. The cases of n = 0.5 are included in Fig. 2. For n = 0.75 and J/t = 3.3, 
3.4 and 3.5, the total energies as a function of variational parameters are shown in Figs. 
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13(a)-(c), respectively. The energy is improved on that of the GWF in every case. The 
value of a giving the energy minimum becomes large with increase of J/t; this means 
the enhancement of the attractive correlation. 

One can read that the aspects of the minima are quite different among three figures 
of Figs. 13(a)-(c). For simplicity, let us consider the case of (3 = 1 (empty circle). In each 
figure the energy curve of j3 = 1 has two local minima, that is, one around a = 1 ~ 2 
(minimum H) and one around a = 15 (minimum S). In Fig. 13(a) the lowest energy is 
given by the minimum H, and the minimum S has a higher value. In Fig. 13(b) the 
two minima have comparable values. On the other hand, for J/t = 3.5 (Fig. 13(c)) 
the situation is opposite to Fig. 13(a). Thus we find a switching from the state of the 
minimum H to one of the minimum S around J/t = 3.4. 

Let us see that the optimal state for J/t > 3.4 represents the phase separation. 
First, the snapshots of electron configurations taken along a Monte Carlo (MC) sweep 
(Fig. 14) show a formation of an electron cluster. Although small density fluctuation 
can be seen near the boundary, this edge effect will be irrelevant as N a — > oo. Second, 
the minimized energy for J/t > 3.4 is quite close to the energy of the Heisenberg chain 
corrected by electron density (Fig. 2). 

To see this more quantitatively, charge- density correlation function in real space, 

jV r = -L^[(n^ +r )-n 2 ], (4.4) 



a 



is plotted in Fig. 15 for f3 = 0.75. For a = (GWF), N r has a negative dip in a short 
distance due to the exchange hole of the original Fermi sea, and falls off as r~ 2 in a 
long distance and converges to zero. When a is switched on, the effect of the exchange 
hole is compensated by the weak attractive correlation; on the other hand, long range 
part of N r hardly changes. As the value of a increases further, N r abruptly changes its 
behavior at a = 10. It approaches the value of the completely separated phase (a = oo) 
expressed as 

N = [n-n 2 -r/N & (0 < r/N a < 1 - n) 1 

iVr \2n-n 2 -l (1 - n < r/N a < 1/2) tOT 2 S S 1 J 

Nr ~ \ -n 2 (n < r/N a < 1/2) < n < ^ (4.5ft) 

In Fig. 16 Ni, N r= N a / 2 and 5/v a /4 = —N a [dN r /dr) for r = A^ a /4 are plotted versus a. 
The transition at a ~ 10 can be seen. When a large cluster of electrons is formed, N± and 
N r=Na / 2 are most influenced by the density fluctuation near the boundary. Inversely 5 r 



14 



for r ~ iV a /4 depends mainly on the fluctuation in the middle of the electron cluster, 
so that it represents the phase transition sharply. 

In Fig. 17 S(k) is shown for some values of a. For a = it has a cusp at k = 2kp = 
ir/2. As a becomes large S(tt) is enhanced, showing antiferromagnetic correlation for 
the Heisenberg chain. For a = 30, S(k) of the AJWF is considerably similar to the 
Heisenberg spin system corrected with n. The behavior of S(k) is also consistent with 
the phase separation. 

In this way we find the phase boundary to the phase separation is at J/t = 3.4 for 
n = 0.75. Similarly we obtain J/t ~ 3.2 (n = 0.25) and J/t ~ 3.3 (n = 0.5). These 
values are in good agreement with the exact-diagonalization result. 

§5. Tomonaga-Luttinger-Liquid-Type Jastrow Function 

In this section we study the trial state proposed by Hellberg and Mele. 23 ) For the 
variational states treated in the previous section, the correlation exponent is always 
the same as the Fermi liquid (K p = 1). This is because the Jastrow factor 7](r) used 
in (4.2) approaches 1 rapidly, namely is short ranged. If the correlation factor t](r) is 
long-ranged, the correlation exponent becomes non-trivial which is consistent with the 
TL liquid behavior. Hellberg and Mele 23 ) showed that a variational state 

\F(r thrjl )\^F, (5.1) 

gives non-trivial exponents, where \F(r^, rji)\ is a Slater determinant of all the electron 
positions. Actually the correlation exponent is related to the variational parameter v 

as 36 ) 

K p = — - — , (5.2) 
p 2z/ + l' v ; 

which was derived using the asymptotic Bethe Ansatz and the scaling relations of the 

TL-liquid theory. 

If we rewrite this wave function in the form (4.1) using the Vandermonde's deter- 
minant identity, the corresponding Jastrow correlation can be written as 



rj(r) 



— sin — r 

7T AL ; 



(5.3) 



When v > 0, this correlation is repulsive and is attractive for v < 0. As we can see from 
(5.3), the factor r](r) behaves as N£ at the longest distance, r = N & /2. It depends on the 
system size Af a , and thus it is a very long-ranged Jastrow factor which is different from 
the conventional Jastrow factor, (4.2). This long-range behavior stands for the effect 
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of phase shift of two-particle collisions and eventually changes the correlation exponent 
from the non-interacting value. 

In Fig. 18, we show the phase diagram obtained using this trial state. The solid 
lines represent the contours of constant correlation exponent K p . Comparing this phase 
diagram with that obtained from the exact diagonalization of small clusters, 11 ) we find 
that the variational wave function is very good in the vicinity of J/t = 2 including the 
correlation exponents. However agreement is not so good near J/t ~ 0. For example, 
the optimized variational parameter v for J/t = is v = 1, (v = 0.75, v = 0.5) near n = 
0, (n = |, n = 1), respectively, while the exact exponent is K p = |, which corresponds to 
v = |, regardless of n. One can see this deviation directly by examining the correlation 
functions. Since the density correlation function N(k) is exactly the same as that of 
the free spinless fermions at J/t = 0, there must not be a peak at 2/c F - However N(k) 
calculated with r\ of (5.1) has a small peak at 2/c F owing to the deviation of the exponent 
(Fig. 19(a)). 

To understand the reason of this deviation, we rewrite the variational state (5.1) 

as 37 ) 



|F(r it ,r i; )r* P , = ni* 


-zj\ v n\z%i 


- Z 3]\\\\Zi[- Zj[\ 


i<j 


i<j 


i<j 


= 115* 


-^rvn 


\Zi] — Zjl\, 


i<j 


T-,3 





apart from the sign. Here Zj a = exp(27rirj a /N a ) and {zj} represents all the electrons 
irrespective of their spins. On the other hand, the exact density correlation function at 
J/t = is reproduced by the free spinless fermion, Yli<j \ z i ~ z j\- We can see that it is 
difficult to reproduce the exact density correlation using the form (5.4). 

In this connection, for the region of J/t ~ 0, a better trial function will be obtained 
as the form, \1/ = x$sf on the analogy of the exact eigenfunction in the small- J/t limit, 
12 ) where $sf is the wave function of spinless fermion and x 1S that for spin degree of 
freedom. 38 - ) 

Next we make a comparison of correlation functions. Figure 19(a) shows n(k), S(k) 
and N(k) for the optimized wave function (5.1) at J/t = and n = 0.5, together with the 
result of exact diagonalization. The reproducibility of the global features is not so good. 
In the same figure we compare the correlation functions calculated in the optimized 
RJWF; the two variational states give almost identical results. Actually we find the 
two optimized forms of r/(r) of two variational states are similar in short distance. This 
suggests that the global features of the correlation functions are determined mainly by 
the short-range behavior of the Jastrow factor. 
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case 2) rj(r) = 



Figures 19(b) and 19(c) show similar comparisons between the variational states 
and the exact diagonalization at J/t = 1 and 2.6, respectively. Since the values of J/t 
is nearer to 2, three results are close to each other except for the vicinity of hp in n(fc). 

Before closing this section we confirm that the long-range part of the correlation 
factor controls the critical behavior of correlation functions, while the short-range part 
determines the global features of correlation functions as well as energy. Let us consider 
the wave functions with four different types of correlation factors, two of which connect 
the correlation factors of two kind of variational states. Namely, 

case 1) rj(r) — 77tll(?") for all r, 

rtFh(r), for r < r c ; 

??TLL(r), for r > r c , 

case3) V{r) = l m ^ ^ r <r c , 

\ VFL{r), for r > r c , 

case 4) rj(r) = ?7fl(^) for all r, 

where ?7tll denotes rj(r) of (5.3), ?7fl either that of the RJWF or of the AJWF (4.2), 
and r c is a certain value of r which divides rj(r) into a short-range part and a long- 
range part. We choose the variational parameters and r c so that the two parts may be 
smoothly connected. By using these correlation factors, we perform VMC calculations 
for the energy, n(k), S(k) and N(k). We fix the system size to N a = 100 for n = 0.5 
and take 5 x 10 4 samples for each case. 

For the repulsive case we choose r c = 5, v = 0.5 (the optimal value for J/t ~ 0.7) 
and ( = 6.9847. In Table III (a) we show expectation values of several physical quantities 
for the four cases. 8n{k) in the last two columns indicates the differences of n(k) between 
the two possible /c-points adjacent to kp from inside and outside, respectively. For the 
bulk quantities like E t /t, Ej/J and n(k = 0), the cases 1) and 3) give similar values, 
which are different from the cases 2) and 4). This implies that the short-range part of the 
correlation factor is responsible for the bulk properties like energy and the magnitude 
of correlation functions. On the other hand, for 8n{k) the cases 1) and 2) have similar 
values compared with the cases 3) and 4). This indicates the long-range part controls 
the critical behavior of correlation functions. 

In Table 111(b), the result of the same quantities in the attractive regime is summa- 
rized. Here we select r c = 10 and v = —0.15 (the optimal value for J/t ~ 2.5), and the 
two parameters are determined by requesting ?7tll(1) = ^aj(I) and ^tll^c) = ^aj(^c) 
as a = 0.6804 and f3 = 0.5481. The tendency agrees with the repulsive case. These 
facts explain the similarity of n(k) and the slight energy difference shown in Figs. 19. 
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§6. Static Properties of Charge and Spin 

In this section we focus on the behavior of charge susceptibility Xc and spin susceptibility 
X s . Since the trial states in the half filling are always Mott-insulating for the t-J model, 
we can see the critical behavior variationally, and discuss the difference between the 
present case and the Brinkman-Rice transition (Appendix B). 

We calculate Xc and Xs from the VMC data for systems with finite sizes from the 
formulae: 

a2 77 t\t2 

Xc 1 =Q^2= -fi E ( N + 2, 0) + E(N - 2, 0) - 2E(N, 0)}, (6.1) 

and 

X7 1 = q^2= ^i E ( N , 8) + E(N, 0) - 2E(N, 4)}, (6.2) 

where E(N,M) is total energy per site of iV-electron systems with M = — iVj_. 
Henceforth we take 1/t as the unit of Xc and Xs- We confirm that the size dependence is 
mostly negligible for the systems we use: N a = 100 ~ 200. In contrast to the Hubbard 
model in the strong coupling regime, relatively accurate estimate of x is possible for 
the ID t-J model because of the less statistical fluctuations. In this section, we mainly 
use the wave function (5.3) and GWF. Typically 5 x 10 4 ~ 2 x 10 5 independent VMC 
samples are used for each value of the parameter v. To search the optimized value of 
E(N, M) for each N and M, we pick out 61~121 values of v between —1 and 2. 

6.1 Charge Susceptibility 

In Fig. 20 Xc is shown for several values of J/t. Here we observe that with increasing 
J/t, Xc becomes large. This is because enhanced attractive interaction between electrons 
enlarges charge compressibility k (= x c /n 2 ). In the high electron density region (n ~ 
1) Xc for every value of J/t is divergent as n — > 1. This divergence is due to the 
strong correlation effect, which suppresses charge fluctuation; this is in contrast with 
the non-interacting case, where x" 1 = 7rsin(n7r/2) and remains a finite value 1/n as n 
approaches 1. For J/t = and 2 the exact values are known; the spinless fermion result 
X^ 1 = 27rsin(n7r) for the former case, and the Bethe Ansatz solution for the latter. 4 ' 9 ) 
In both cases the exact Xc diverges as Xc oc 1/6 (5 = 1 — n) in the limit of n — > 1. The 
variational result is quantitatively similar to the exact one. To see this divergence more 
closely, we plot chemical potential = dE/dn vs. 5 2 for J/t = and 2 in Fig. 21. For 
J/t < 2 we can fit \x as \i = fio — a ^ 2 m the vicinity of the half filling, although the 
linearity is not so clear for J/t > 2. Since Xc = dn/d/j,, the coefficient a is related to 
the divergence of Xc as Xc = 2a/5(5 — > 0). From Fig. 21, we determine this coefficient 
2a for a couple of values of J/t, which are shown in Table IV. 
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As can be seen in Table IV, the estimation by the variational state is in good 
agreement with the exact value at J/t = 0. However at J/t = 2, the deviation is quite 
large. This is probably because the range of linear behavior of fi vs. 5 2 is rapidly 
reduced as J/t increases. 9 ) 

The divergence of Xc near the Mott transition has been understood from the Bethe 
Ansatz. Using the correlation exponent K p and the charge velocity v c , Xc can be written 
as 5_8) 

2K P 
Xc = -■ 

TV V c 

Since v c vanishes linearly as 5 — > 0, v c = b8, Xc diverges as 39 ' 10 ' 14 -' 

tv b 5 

Since K p — > 1/2 as 5 — > 11 ) the coefficient 2a is solely related to the coefficient b, which 
is 2tt > b > 3yrC(3)/(16(ln2) 2 ) for < J/t < 2. 11 ' 9 ) However, for J/t > 2, we cannot 
expand fi as fx = fio — ad 2 , but fit fi as fi = fi$ — ad p instead. The leading power p 
is estimated as 1 < p < 2, namely p = 1.8 (for J/t = 2.0), 1.5 (2.5) and 1.3 (3.0). 
Although Xc diverges as 5 — > for J/t > 2, we have not understood the origin of this 
power. This may be due to the poorness of the trial wave function of (5.1) or due to 
an anomalous dependence of v c as v c ~ a5 p_1 . At any rate, we can consider that this 
divergence is attributed to the divergence of density of state at the band edge of the 
spinless fermion. 

In low electron density area (n ~ 0) Xc again diverges except for the case of J/t = 
2.5, in which the variational state becomes unstable against phase separation and Xc 
becomes negative. One can find that the divergence of Xc is proportional to 1/n by 
plotting fi as & function of n 2 . The coefficients of 1/n are shown also in Table IV. Since 
for J/t = 2 the trial wave function (GWF) is exact in the limit n — > 0, the exact value is 
within the range of error. Furthermore, the GWF result and the exact one are extremely 
close to the non-interacting gas, which means a "free electron" state is realized in the 
supersymmetric case. Actually using the variational energy (3.12), we obtain 

7T 2 n 

which is the same as in the non-interacting case. As n increases, however, the sup- 
pression of charge fluctuation in the t-J model prevents the wave function from the 
free-electron behavior. 

For J/t = the variational value is a little different from the exact one. This 
difference corresponds to the deviation of K p in this described in the previous 
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section. The divergence of Xc is not due to the correlation effect, but simply to the 
divergence of the density of state at the band edge. 

6.2 Magnetic Properties 

First, we discuss spin susceptibility under zero field. Figure 22 shows Xs for a couple of 
values of J/t. In sharp contrast to Xo Xs does not diverge as n — > 1 for every value of 
J ft but remains a finite value close to that of the Heisenberg chain. 40 '' The divergence 
of Xc hardly affects % s ; this is due to the separation between spin and charge degrees of 
freedom in the low-energy excitations. 

For J/t = 0, all the spin configurations are degenerate, hence Xs is infinite. This 
aspect is special in ID. By introducing J/t, this degeneracy is lifted and the value of 
X s becomes finite except for the low density limit, which is shown in Fig. 22. In the TL 
liquid theory, Xs can be represented as 5 ' 6 ' 39 ' 7 -' 

2 1 

Xs 

7T V s 

where v s is the spin velocity. Basically v s is proportional to the exchange coupling J and 
thus Xs decreases as J/t increases, which means that the enhanced exchange coupling 
hinders the response of spins to the magnetic field. For a low density area (n ~ 0) Xs 
diverges since v s oc n, which is due to the divergence of density of states. On the other 
hand as n — > 1, Xs approaches that of Heisenberg chain, 

2 1 

Xs -> 

TT fsw 

with v sw being the spin wave velocity of Heisenberg chain: v sw = J7t/2. 

For the supersymmetric case, we observe that the variational results agree quite 
well with that of the non-interacting system (Pauli paramagnetism) x^ 1 = 7r sin(n7r/2), 
especially in low density region. Here the idea of the GWF as a "free electron" state 
is again useful. Like Xc however, the two results become a little apart as n — > 1. For 
J/t = 2.5 and n < 0.1, the decrease of Xs is due to the phase separation. 

Now let us turn to the case under a finite magnetic field gn^H, and here we put 
g = 2 and ^3 = 1. Applying an external field along z-axis, a Zeeman term, 

H ext = -2H Sj = ~N a Hm, (6.3) 

3 

is added to the original Hamiltonian (2.1). Total energy per site is written as, 

E = E(m,n) - Hm, (6.4) 
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where E(m, n) is total energy per site with electron density n and magnetization m 
under zero field. 

In Fig. 23 E(m, n) is shown for n = 0.5 and 1.0 as a function of m (0 < m < n). 
For each case the energy is a monotonically increasing function of m. By minimizing 
the total energy (6.4), we obtain the magnetization curves, which are shown in Fig. 24. 
The critical field H s at which spin saturates is determined from the slope of E(m,n) 
near m = n. The value of H s is 0.34(J/t = 1.0), 0.99(2.0) and 1.32(2.5) for the quarter 
filling, and 2.000(2.0) for the half filling. 

For the half filling, the data of the GWF is very close to the exact value (Heisenberg 
antiferromagnet) 40 ) for all the range of if. 41 ) In the weak field limit, m/H is nothing 
but x s , the similarity of the two results is obvious from Fig. 22. 

For the quarter filling, m saturates at smaller H as J/t decreases. This is naturally 
understood since the energy to excite the system to higher spin state becomes less as 
J/t decreases. In the supersymmetric case, the GWF is in good agreement with the BA 
result also for the quarter filling. 42 ) Furthermore, the non-interacting result agrees well 
for all the values of H/t. 

6.3 Comparisons and Discussions 

First, we compare the above result with the Hubbard model. Exact results 43 ' 44 ' 39 -* show 
that in the limit of n — > 1, Xc diverges as Xc = ot/5, where a is a numeral factor which 
depends on U/t. 39 ) The value of a changes from zero for U/t = to 1/2-7T 2 for U/t = oo, 
which is the same with the t-J model with J/t = 0. On the other hand, Xs converges 
upon a finite value as n approaches l. 43 -* In the strong-t/ limit the expansion coefficient 
in t 2 /U of x s is the same as the t-J model. These features qualitatively agree with the 
variational results of the present study as well as the exact one for the t-J model. 

Meanwhile, Otsuka 45 -* and Furukawa and Imada 10 -* investigated the critical behavior 
of Xc and Xs for the Hubbard model on 2D square lattice by using quantum Monte Carlo 
methods. According to their results, essentially the same properties with the ID models 
are observed. Thus one can consider that a universal profile of a kind of Mott transition 
appears in these results of Xc and Xs- 

Keeping these in mind, next let us compare with the Brinkman-Rice transition 
in the Hubbard model. As summarized in Appendix B, Xs calculated with Gutzwiller 
approximation diverges when U — > U c (n = 1) or n — > 1 (U > U c ). On the other hand, 
Xc remains finite for a finite value of U, even if U > U c , although effective mass m*/m 
diverges similarly as Xs- This indicates that the Brinkman-Rice transition describes a 
quite different type of metal-insulator transition from those in the ID and 2D Hubbard 
models and the ID t-J model. 



21 



Lastly we mention the results in the long-range t-J model. 30 ) In this model, sus- 
ceptibilities are written as 



Xc _1 = 



X7' = y(l-m). 



(6.6) 



Notice that Xs does not depend on n and Xc does not on m, since the contributions 
to energy of n and m are mutually independent, and that there are no system size 
dependence, because the size-dependent terms in energy are exhausted in the linear 
order. The value of Xs is constant (2/-7T 2 ) under zero field irrespective of n, which is the 
same for the nearest-neighbor Heisenberg model and is close to the variational value 
(GWF). The divergent behavior of Xc near the half filling is also similar to the t-J 
model (see Table IV). On the other hand, a quite different feature appear in the low 
electron density region, where there is no divergence. This difference is originated in 
the band structure of the model; the non-interacting long-range model has linear energy 
dispersion and there is no divergence in density of state at the band edge. From this 
we can see that electron correlation affects severely high-electron-density regime, while 
in the low-density regime the density of states of the original non-interacting system 
determines the charge susceptibility. 

In Fig. 24 we also plot the magnetization curve of the long-range t-J model: 46 ' 32 ) 



This formula does not depend on n. The critical field H s for the long-range model is 
somewhat larger (7r 2 /4 for n = 1 and 37r 2 /16 for n = 0.5) than the value of the ordinary 
t-J model. And the difference is larger for n = 0.5. This is probably because the 
long-range exchange terms tend to disturb ferromagnetic spin alignment. Especially in 
the low density region, where the particle distance is large, the long-range terms play 
important roles. 




(6.7) 
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§7. Summary and Discussion 

We have pursued the ground-state properties of the one-dimensional t-J model in the 
light of wave functions. By comparing the variational Monte Carlo results with those 
of the exact diagonalization, the Bethe Ansatz and the Gutzwiller approximation, we 
have obtained some remarkable aspects as follows: 

[1] From diagonalization, exact energy, momentum distribution, and spin and charge 

correlation functions are obtained, which show unusual behaviors as J/t increases. 

[2] In the supersymmetric case (J/t = 2) the Gutzwiller wave function is an extremely 

good state for bulk quantities. In the low-density limit the GWF becomes exact. 

[3] Whole parameter space spanned by electron density n and coupling strength J/t 

is well described by the repulsive or the attractive Jastrow-type wave functions. The 

parameter space can be roughly divided into: 

(a) J/t < 2 repulsive region (RJWF), 

(b) J/t = 2 free-like region (GWF), 

(c) J/t > 2 attractive region (AJWF). 

In the region (a), the electron hopping term is dominant. In the region (c), dominant 
is the exchange term, which induces attractive interaction. In the supersymmetric case 
(b), the two terms are well balanced and a kind of "non-interacting" state is realized, 
especially for the low electron density. 

[4] Phase transition from a homogeneous state to a phase-separated state is quantita- 
tively described within the AJWF. 

[5] In the Jastrow wave functions, a short-range part of intersite correlation factors are 
responsible for relatively high energy processes which determine bulk properties like 
energy and magnitude of correlation functions. On the other hand, a long-range part 
is mainly concerned in low energy processes near the Fermi surface, which cause the 
critical properties characteristic of the Tomonaga-Luttinger liquid. 
[6] The Jastrow wave functions reproduces charge and spin susceptibilities and magne- 
tization curve correctly, in contrast with the Gutzwiller approximation. 

Keeping these results in mind, we mention some remaining issues. 

In the region of low electron density and J/t > 2, there exists a spin gap state 
(without a charge gap). 11 ' 47 -' We have not found an indication of a spin gap in the trial 
functions used. On the other hand, Chen and Lee 48 ) introduced a trial state for a gas 
of singlet pairs, and showed that there is a region where this function is stabler than 
the TL-liquid wave function (5.1). 

An interesting extension of the present method is to 2D systems. In 2D we do not 
know even the ground state — the Fermi liquid or the TL liquid for the metallic regime. 
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Furthermore, magnetically ordered phases can be stabilized, near the half filling. 49 ) 
Actually the TL-liquid wave function has been extended to a 2D system by Valenti and 
Gros. 50 ) According to their results, the energy lowering by their function is very small 
(1%) compared with the simple Gutzwiller wave function. Also the critical exponent 
at k-p in the momentum distribution is small, assuming that it exists. Since critical 
properties are in a low energy scale, it may not be easy to judge the realization of a 
TL-liquid state in 2D only by the stability in energy. 

Recently the behaviors of specific-heat coefficient, effective mass, Xc Xs, etc. are in- 
vestigated experimentally for the high-T c superconductors and related Mott insulators, 51 ^ 
in connection with the metal-insulator transition. The results of these experiments to- 
gether with the theoretical calculations 45 ' 10 ^ have suggested reconsideration of the ap- 
propriate Hamiltonian, namely whether the simple t-J model is pertinent to describe 
various aspects of the high-T c cuprates. 

In these contexts, 2D systems have to be further studied with the variation theory. 
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Appendix A: Analytical Approach to the Gutzwiller Wave Function. 



In this appendix, we summarize the analytic expressions for various expectation values 
in the GWF. They were developed by Metzner and Vollhardt, 18 -* and by Gebhard and 
Vollhardt 19 ) for the Hubbard model. In order to apply to the ID t-J model, we carry 
out the calculation of (niUj) and (SfSj). 

The momentum distribution for spin a is obtained by an infinite summation 18 ) 

n a (k) =n° a (k)-(l-g) 2 n_ a n a (k) 



+ 



(1 + 9Y 



£ (9 2 - l) m {l " (1 " 9 2 K(k)}fmAk), (A.l) 



m=2 



where g is identical to i](0) in our notation (eq. (4.1)), and thus g = for the t-J 
model. n^(k) is the momentum distribution of the non-interacting system given by 
n®(k) = 9(kp — \k\), and n a = n/2. In ID case, f m (k) is given by a polynomial of order 
< to. According to the notation in Ref. 18 ), we can summarize as follows: 



( n m R m {k) , 

(in region (I) : < k < kp) 



fm(k) = < 



nrn Qm(k) + C m — 1 , 

(in region (II) : k-p < k < min(3/cF, 27r — 3/cf) 



(A.2) 



a 



m—l 



(in region (III) : 3/cp < k < 7r,when 3/cp < tt) 



n m {Q m {k) + Qm(2ir - k)) + C m _i , 

(in region (IV) : 2n — 3/cf <k<ix, when 2n — 3/cf < tt) 



where 



C m -1 = (-1)' 



2m' 



m d(j') 



^ O ij) k 

«-(*> = E^ss 



4^ ' 



(A.3) 



The coefficients of these Taylor series are determined via a recursion relation; 

2 ( m _ j + i)^) +1 = _(2to - 2j + 1)R$ - Q^XW 
Q^+f = 2(to - 2 J )Q^1 ) + 2mR$ - 4j(m - j + l)(i#> + i^ } +1 ), 



25 



„( m +D fO m + 1 = odd . . 

^ = {-2Q%?\ m + l = even ( ) 

with initial values 

d(j)_ X x 
- ~2 i' ' 

Apparently the series in (A.l) gives an expansion with respect to n in the low 
density region. For m = 2, we get 

1.2 



( ^( k 

2/1. qx 2 „2 



/2(fe) 



(A.5) 



-f(i-I) +f - (II) 

£ in (III) 

V 2 12™ 4 ) 2 I 2™ 4 J + 4 m ^1 V j 

Then the kinetic energy in the dilute limit is 

E t = -2t(n-^- ^n 3 - + 0(n% (A.6) 

Compact analytic expressions for S(k) and N(k) are given in (3.17) and (3.18). 
Their Fourier transforms give the spin and charge correlation functions in real space. 
Through a straightforward calculation, we get 

A{SlS z l+r ) = ^{Si(Trr) - Si((l - n)yrr)}, (A.7) 

and 

(nini +r ) = n 2 + * 9 (cos4fc F r - l) 

, ^ (A-8) 
- — { h (1 - n) cos n7rr) { Si(yrr) - Si( (1 - n)7rr) } , 

where r is an integer, r > 1. When r = 1, (A.7) and (A. 8) give the exchange energy 
given in (3.3). In the low density limit, we obtain 

n 2 n 3 . 4 



(^ +1 ) = - Y -- + 0(n 4 ), 



and 



Total energy is 



9 S 

n n 



(n l n l+1 ) = — + — + 0(n*). (A.9) 



£■ = -2tn +{t- J -)n 2 + ^-tn 3 + (\t -\j)n 3 + 0(n 4 ). (A.10) 



At J/t = 2, we get (3.12). 
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Appendix B: % s and Xc around Brinkman-Rice Transition. 

Brinkman and Rice discussed the metal-insulator transition for the Hubbard model 
based on the Gutzwiller approximation (GA). 52 -* Although it was confirmed by solving 
the variation problem accurately that this transition does not exist in the realistic (1,2 
and 3) dimensions, 17 > 18 ) the conception of the Brinkman-Rice transition is still widely 
used. In this Appendix we briefly review the behavior of Xs and Xc in the GA, especially 
around the Brinkman-Rice transition for the comparison in §6. 53 - ) 

According to the GA, the state of n 7^ 1 in the Hubbard model is always metallic 
and has a Fermi surface with a finite discontinuity q of n(k) at k = kp. On the other 
hand, in the half filling the GA gives a metallic state for U < U c and an insulating 
one for U > U c , where U c = 8\sq\ and £0 is the energy of the non-interacting system; 
the Brinkman-Rice transition occurs at U = U c . When U increases from under U c 
fixing n = 1, Xs as weu as effective mass m*/m(oc g -1 ) diverges as 1/(1 — (U/U c ) 2 ). 
Xs and m*/m remain infinite for U > U c . Meanwhile, Xc decreases with increasing 
U and vanish at U = U c and then increases for U > U c . In Fig. 25 we actually plot 
the numerical GA values of Xs and Xc- According to the Fermi liquid theory, charge 
susceptibility is related to m*/m as 

m*/m Q 

Xc - 1 + F s Xo (tf- 1 ) 

where Fq is the usual Landau parameter and Xc * s the value for the non-interacting 
case. Therefore Fq is more divergent than m* jm at U = U c . 
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Table I Ground-state energies of the ID t-J model for the quarter-filled case (n = 1/2). 

They are extrapolated to the thermodynamic limit using the formula (2.6). and 
e'^ are obtained from the two set of boundary conditions (see the text). The unit 
of the energy is t. 



11+ 

J 1 1 


p a' 

c oo too 


u.u 


— U.DoDDiyi — U.DoDDiyo 


u.z 


— U.ODfofOU — U.0D(DU4o 


n a 
U.4 


— U.0oU4Uzy — U.DoUDZDy 


U.D 


n 70/117R7 n 70/117R8 


U.o 


n 7ooi 7ni n 7001 7ni 


l.U 


— V.tOOOOot — U./OOoOoD 


1 9 
± .Z 


n 789791 1 n 7897900 

— u./oz/zii — u. loz/zuy 


1.4 


-0.8112399 -0.8112396 


1.6 


-0.8409031 -0.8409028 


1.8 


-0.8717049 -0.8717046 


2.0 


-0.9036477 -0.9036477 


2.2 


-0.9367449 -0.9367452 


2.4 


-0.9710250 -0.9710259 


2.6 


-1.0065398 -1.0065415 


2.8 


-1.0433801 -1.0433827 


3.0 


-1.0817130 -1.0817179 


3.2 


-1.1219079 -1.1219533 
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Table II Ground-state energies of the 
quarter-filled case in Table I. 



ID t-J model for n = 3/4 similarly obtained as in the 



J/t 


e e' 


0.0 


-0.4501050 -0.4502187 


0.2 


-0.5216656 -0.5219004 


0.4 


-0.5943201 -0.5944475 


0.6 


-0.6678961 -0.6680324 


0.8 


-0.7422736 -0.7424195 


1.0 


-0.8173650 -0.8175213 


1.2 


-0.8931052 -0.8932726 


1.4 


-0.9694457 -0.9696252 


1.6 


-1.0463508 -1.0465437 


1.8 


-1.1237962 -1.1240038 


2.0 


-1.2017671 -1.2019909 


2.2 


-1.2802580 -1.2805002 


2.4 


-1.3592732 -1.3595368 


2.6 


-1.4388286 -1.4391186 


2.8 


-1.5189551 -1.5192801 


3.0 


-1.5997064 -1.6000835 


3.2 


-1.6811763 -1.6816414 


3.4 


-1.7635421 -1.7641837 
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Table III Comparison of expectation values for four kinds of correlation factors hybridized 
(a) between the RJWF eq.(4.2) (C = 6.9847) and eq.(5.1) {v = 0.5) for the repulsive 
case, and (b) between the AJWF eq.(4.2) (a = 0.6804, /3 = 0.5481) and eq.(5.1) 
{y = —0.15) for the attractive case. Digits include some statistical errors. 



(a) 


Et/t 


Ej/J 


n(k = 0) 


Sn(k) < kp 


Sn(k) > kp 


case 


1) 


-0.60865 


-0.13972 


0.73826 


-0.02925 


-0.02948 


case 


2) 


-0.60959 


-0.12047 


0.70705 


-0.02665 


-0.02811 


case 


3) 


-0.60880 


-0.13955 


0.73804 


-0.01958 


-0.01989 


case 


4) 


-0.60983 


-0.12037 


0.70816 


-0.02242 


-0.02178 



(b) 


Et/t 


Ej/J 


n(k = 0) 


Sn(k) < k-p 


Sn(k) > kp 


case 1) 


-0.55403 


-0.17331 


0.74876 


-0.00589 


-0.01103 


case 2) 


-0.54695 


-0.17619 


0.74777 


-0.00309 


-0.00559 


case 3) 


-0.55392 


-0.17334 


0.74871 


-0.00096 


-0.00553 


case 4) 


-0.54701 


-0.17615 


0.74745 


-0.00090 


-0.00279 



Table IV Coefficient of 1/ra and 1/5 of charge susceptibility in the limit of n — > and 5^0, 
respectively. For comparison we also show the values of the free electron system 
and the supersymmetric t-J model with long-range exchange and transfer, denoted 
by 'Free' and 'L.-r. t-T respectively in the first column. 





n -> 


6^0 


J/t 


Variational 


Exact 


Variational 


Exact 




1 

2 


0.046 
0.053 

0.204 


+ = 0.05066 
A = 0.20264 

7T 


0.054 
0.062 

0.070 


+ = 0.05066 
= 0-21598 


Free 




\ = 0.20264 




Finite Xc 


L.-r. t-J 


Exact 


Finite Xc 


Exact 


4 = 0.20264 



33 



Figure Captions 

Fig. 1 (a) The momentum distribution function, (b) spin correlation function, and (c) 
charge correlation function for the quarter-filled case at J/t = 0.5, 1.0, 2.0 and 3.0 
obtained in the exact diagonalization of the systems with 4, 8, 12 and 16 sites. For 
J = 0, we show the result in the wave function for the U — > oo Hubbard model 
with 52 sites. 12 ' 15) 

Fig. 2 Ground-state energy of the ID t-J model as a function of J/t at n = 0.5. Extrap- 
olated values (N a — > oo) of the exact diagonalization, the energy of the fully phase 
separated (Heisenberg) state, and variational energies of the three types of varia- 
tional functions are compared. The transition to a phase separated state occurs at 
J/t ~ 3.3 shown by a hollow arrow. In the VMC calculations systems with 50-100 
sites are used. 

Fig. 3 Comparison of the total energy per site between the Gutzwiller wave function 
(closed circles) and Bethe Ansatz (Solid line) for J/t = 2 as a function of electron 
density. Although the VMC results for N a = 102 is plotted in the figure, the 
analytic expression for the Gutzwiller wave function is available as in the text. The 
non-interacting case is also shown. 

Fig. 4 Comparison of momentum distribution function between the GWF (solid line) and 
the exact diagonalization (open and closed circles) at J/t = 2. The analytic expres- 
sions for the GWF have been obtained in Ref. 18 ' 19 -* and given in Appendix A. Here 
we show the results of VMC calculations in 60- and 72-site systems with 5 x 10 4 
samples. For the diagonalization, the data in N a = 4, 8, 12 and 16 sites for n = 0.5 
and N a = 8 and 16 sites for n = 0.75 are shown. 

Fig. 5 Comparison of spin and charge correlation functions between the GWF (solid line) 
and the exact diagonalization (open and closed circles) at J/t = 2. Numerical 
results are obtained similarly as in Fig. 4. 

Fig. 6 Correlation factor r](r) for three types of Fermi-liquid-type wave functions. Typical 
values of the parameters are chosen for each case. 

Fig. 7 Variational expectation values of total energy in the supersymmetric case for (a) 
the RJWF and for (b) the AJWF. The value of the GWF is shown by an arrow on 
the vertical axis. In this calculation 3 x 10 4 samples are used for a 60-site system. 

Fig. 8 Expectation values of two energy component Tit (solid circle) and Tij (open circle) 
per site for the RJWF for n = 0.5. The value of the GWF (£ = 0) is shown by 
arrows on the vertical axis. We use 3 x 10 samples for a 60-site system. 

Fig. 9 Energy expectation values for the RJWF at J/t = 1.0. Arrow on the vertical axis 
is the value of the GWF and that of the exact diagonalization (extrapolated values 
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for iV a — > oo). 

Fig. 10 Expectation values of two energy component H t (solid circle) and H j (open circle) 
per site for the RJWF and for small-C region. In this VMC calculation 10 5 samples 
are collected for a 72-site system. 

Fig. 11 Discontinuity of momentum distribution n(k) at hp as a function of n. The op- 
timized Fermi-liquid- type wave function is used for each value of n and J/t. The 
size of the symbols represent the relative magnitude of possible error. For the 
cases n — >■ and J/t — > and for the AJWF it is not easy to determine the value 
accurately owing to the difficulty in optimizing the parameters. 

Fig. 12 Variational expectation values of (a) H t and (b) Hj per site for the AJWF. The 
arrow on the left vertical axis shows the value of the GWF (a = 0). The arrows 
on the right axis in (b) shows the value of the Heisenberg chain (Bethe Ansatz) 
multiplied by n. The sample number is 3 x 10 4 for each point and N a = 72. 

Fig. 13 Total energy of the AJWF for a couple of values of J/t around the phase transi- 
tion. The same symbol with Fig. 12 is used for each value of f3. (a) corresponds 
the homogeneous regime, (b) is the case near the phase transition, and (c) is in 
the phase-separated regime. The arrow on the left vertical axis shows the value 
of the GWF. The Heisenberg value corrected by the electron density is given by 
-n(J/t) In 2: -1.8195 • • • for n = 0.75 and J/t = 3.5. 

Fig. 14 Snapshots of the electron configurations in the VMC sweep for N a = 60, n = 0.5 at 
J/t = 3.4 (separated phase). Each horizontal line represents the one-dimensional 
system. An solid (open) circle means an up- (down-) spin electron, and an null 
space (horizontal line) an empty site. In the Monte Carlo sweep, the configuration 
of the system evolves vertically. Before taking these snapshots, 3000 MC steps are 
discarded for obtaining the statistical equilibrium. Sampling interval is determined 
(/ = 20MCS) so as to make the acceptance ratio per electron larger than the unit. 

Fig. 15 Real-space charge density correlation function for the AJWF are shown for various 
values of a. The value of (3 is fixed at 0.75, which gives minima approximately 
around the phase transition. The solid line shows the values of the completely 
separated phase (a = oo). The system size is 60 sites and 5 x 10 samples are used. 

Fig. 16 Real-space correlation function N r for the nearest neighbor sites and its derivative 
5 r for r = N a /2 are plotted as a function of a. Arrows on the vertical axes of 
the left and right represent the values of a = (GWF) and a = oo (completely 
separated state), respectively. The used data are the same ones in Fig. 15. 

Fig. 17 Spin correlation function for the AJWF for various values of a. The solid line 
represents the value of the Heisenberg chain obtained by the exact diagonalizations, 
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which is corrected by the electron density. The system size is 60 sites and 5 x 10 4 
samples are used. 

Fig. 18 Phase diagram of the ID t-J model calculated in the TL-liquid wave function 
(5.1). The curves show the contours of constant correlation exponent K p . The 
used system has 100 sites. 

Fig. 19 Comparison of momentum distribution function n(k), spin and charge correlation 
functions, S(k) and N(k) among the TL-liquid wave function (bold solid line), the 
Fermi-liquid-type wave function (dot) and the exact diagonalization (open diamond, 
open circle and exceptionally thin solid line for N(k) of J/t = 0) for (a) J/t = 
(u = 0.75 for the TL-liquid state, C = 1.6 for the RJWF), (b) J/t = 1 (u = 0.37, 
C = 0.7) and (c) J/t = 2.6 (u = -0.19 for the TL-liquid state, a = 0.4 and 
(3 = 1.0 for the AJWF). For the variational calculations, 10 5 samples are averaged 
for 220-site systems. Total energy for each case is also given in digits (the last 
digit for each include the ambiguity due to the statistical fluctuations). For exact 
diagonalization, 52-site system is used in (a) and up to 16-site systems in (b) and 
(c). 

Fig. 20 Charge susceptibility vs. n for some values of J/t. Symbols are the results of the 
optimized TL-liquid state. Solid line for J/t = and 2 9 ) represents the exact 
analytic value. Dashed line is the result for the non-interacting system. Dotted 
line for J/t = 1.0 and 2.5 is a guide to the eyes. The sizes of the symbols represent 
the relative magnitude of possible error. 50 ~ 210-site systems are used. 

Fig. 21 Chemical potential as a function of 5 2 for J/t = and 2. Symbols indicate different 
system sizes, namely N a = 200 (solid diamond), 150 (open circle), 100 (solid circle), 
70 (solid square), 50 (open diamond). 

Fig. 22 Spin susceptibility vs. n for some values of J/t. Open circle is the result of the 
GWF for J/t = 2 and open diamond is the result of the optimized TL-liquid state. 
Solid line for J/t = 2 represent the exact analytic value. 4 ' 9 -* Dashed line is the value 
for the non-interacting system. Dotted line for J/t = 1.0 and 2.5 is a guide to the 
eyes. Arrow on the right vertical axis shows the exact value for the Heisenberg 
antiferromagnet 2/n 2 for J/t = 2. 40 ) The sizes of the symbols represent the relative 
magnitude of possible error. 50 ~ 210-site systems are used. 

Fig. 23 Total energy under zero field as a function of m calculated with the TL-liquid wave 
function. For the half filling, where the TL-liquid state is reduced to the GWF, 
the value in the supersymmetric case is shown for an example. For n = 1.0 (0.5), 
a 102-(100-)site system is used. 5 x 10 samples are averaged. 

Fig. 24 Comparison of magnetization curves for n = 0.5 and 1.0. Symbols represent the 
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result of the variational functions. The optimized TL-liquid wave functions are 
used except for the supersymmetric case, in which the GWF is substituted. Solid 
line is the Bethe Ansatz result for J ft = 2. 40 ' 42 ^ Dotted line is the result for the 
non-interacting system. The result for the supersymmetric long-range t-J model is 
shown by dashed line. 46 ' 32 ) 
Fig. 25 (a) Charge and (b) spin susceptibilities calculated with Gutzwiller Approximation 
for the Hubbard model as a function of electron density for some values of U/U c . 
Inset in (a) represents Xc vs. U/U c for the half filling. The ID cosine band is 
assumed. 
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